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Geometric phases in electric dipole searches with trapped spin-1/2 particles in general 
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The important role of geometric phases in searches for a permanent electric dipole moment 
of the neutron, using Ramsey separated oscillatory field nuclear magnetic resonance, was first 
noted by Commins [Am. J. Phys. 59, 1077 (1991)] and investigated in detail by Pendlebury et al. 

[Phys. Rev. A 70, 032102 (2004)]. Their analysis was based on the Bloch equations. In subsequent 
work using the spin density matrix Lamoreaux and Golub [Phys. Rev. A 71, 032104 (2005)[ showed 
the relation between the frequency shifts and the correlation functions of the fields seen by trapped 
particles in general fields (Redfield theory). More recently we presented a solution of the Schrodinger 
equation for spin-1/2 particles in circular cylindrical traps with smooth walls and exposed to ar¬ 
bitrary fields [Steyerl et al, Phys. Rev. A 89, 052129 (2014)]. Here we extend this work to show 
how the Redfield theory follows directly from the Schrodinger equation solution. This serves to 
highlight the conditions of validity of the Redfield theory, a subject of considerable discussion in 
the literature [e. g., Nicholas et. al., Progr. Nucl. Magn. Res. Spectr. 57, 111 (2010)]. Our results 
can be applied where the Redfield result no longer holds, such as observation times on the order 
of or shorter than the correlation time and non-stochastic systems and thus we can illustrate the 
transient spin dynamics, i.e. the gradual development of the shift with increasing time subsequent 
to the start of the free precession. 

We consider systems with rough, diffuse reflecting walls, cylindrical trap geometry with arbitrary 
cross section, and field perturbations that do not, in the frame of the moving particles, average 
to zero in time. We show by direct, detailed, calculation the agreement of the results from the 
Schrodinger equation with the Redfield theory for the cases of a rectangular cell with specular walls 
and of a circular cell with diffuse reflecting walls. 

PACS numbers: 28.20.-v, 14.20.Dh, 21.10.Tg 

I. INTRODUCTION 

The experimental search for a permanent electric 
dipole moment (EDM) of the neutron or other particles is 
motivated by its potential to search for deviations from 
the Standard Model. The most sensitive experimental 
approach so far, for neutrons, uses ultracold neutrons 
(UGNs) trapped in a measurement cell together with co¬ 
magnetometer atoms like ^®®Hg. They are exposed to a 
weak static magnetic field Bq and a strong electric field 
E applied parallel or anti-parallel to Bq. Employing the 
Ramsey method of separated oscillatory field magnetic 
resonance the experiment is designed to search for small 
changes of the Larmor frequency due to the presence of 
the electric held. 

At the high level of sensitivity already achieved in pre¬ 
vious EDM experiments m and expected to be im¬ 
proved in current EDM projects (e. g., m\) a signifl- 
cant source of potential systematic error appears to be 
geometric phases. Their role in EDM experiments was 
flrst noted by Gommins [1] and analyzed by Pendlebury 
et al. [3] on the basis of the Bloch equation for spin evo- 
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lution in a time-dependent magnetic held. Subsequent 
work by Lamoreaux and Golub m and Barabanov et al. 
[a used the Redfleld approach [T2UI4] based on the spin- 
density matrix and extended the scope of such studies to 
general NMR physics. The original magnetic held model 
with a constant vertical gradient |3] of the Larmor held 
was extended in Refs. pAUIS] to include more general 
macroscopic flelds and the microscopic held of a mag¬ 
netic dipole. 

In Ref. [19] we showed that the frequency shifts due 
to geometric phases can also be determined by directly 
solving the Schrodinger equation to second order in the 
perturbation. Solutions were presented for circular cylin¬ 
drical traps with smooth walls and arbitrary macroscopic 
magnetic fields as well as for the microscopic field of a 
magnetic dipole. 

The direct solution of the Schrodinger equation in¬ 
cludes vertical spin oscillations and the non-stationary, 
transient spin dynamics, i.e., the gradual development of 
the shifts for an arbitrary number of wall reflections sub¬ 
sequent to the start of free spin precession in the Ramsey 
scheme. 

In the present article we extend this approach to in¬ 
clude diffuse reflectivity from rough walls and to vertical 
cylindrical measurement cells with arbitrary cross sec¬ 
tion. Their top and bottom are assumed to be perfectly 
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flat. In Sec. we develop an analytic solution of the 
Schrodinger equation for perfectly diffuse reflection, in 
the xy plane only, on the cylindrical walls, both for a 
constant gradient field in circular cylindrical geometry 
and for a more general field model where Bg and E are 
not exactly aligned. 

Our work is the first to address directly the effects of 


diffuse wall reflections in different geometries. In Sec. Ill 


we consider a generic cross section of the measurement 
cell and calculate the frequency shifts (that linear in E 
and those quadratic in electric field or magnetic field) for 
cells in the form of a rectangle with arbitrary aspect ratio. 
Rectangular cells are foreseen for EDM project [7]. Our 
results are valid over the full range of in-plane particle 
velocity v for given Larmor frequency wg. We compare 
them with the simulations described in Sec. IIVI and ob¬ 
serve consistency between the two approaches. The non- 
adiabatic and adiabatic limits of frequency shift agree 
with the universal predictions of Refs. [HUHHH]. In 
Sec. 1^ we establish the equivalence of our method with 
other approaches nnmiiin]- The agreement is notable 
in view of the different sets of assumptions made by each. 

The usual approach to these problems, the Redfield 
theory of frequency shifts and relaxation based on the 
von Neumann equation for the density matrix, makes 
use of a plethora of assumptions [12] and has resulted in 
a considerable literature concerning its range of validity, 
e. g. [22]. The density matrix represents the effects of 
both the quantum mechanics and statistical fluctuations 
of the system under study. In our approach, which in¬ 
volves first solving the Schrodinger equation, to second 
order in the perturbation, for the wave function for an 
arbitrary time varying field and then applying the statis¬ 
tical averaging to the result, there is a clear separation 
between the quantum and statistical effects. In fact our 
solution can be easily applied to cases where the Redfield 
result no longer holds, such as observation times shorter 
than the correlation time, and non-stochastic systems. 


II. FREQUENCY SHIFT FOR A CIRCULAR 

CYLINDRICAL TRAP WITH ROUGH WALL 

A. Fundamentals 

In the EDM experiments the cell walls are never per¬ 
fect mirror surfaces, thus over the course of some 10^ to 
10^ consecutive reflections randomization of trajectories 
is established for any initial distribution. The total num¬ 
ber of wall reflections occurring during the measurement 
period in EDM experiments is much larger. 

In the present article we neglect scattering on gas 
molecules and consider perfectly diffuse surface scatter¬ 
ing in the xy plane perpendicular to the vertical sym¬ 
metry axis of the cell (Fig. [^, thus it suffices to analyze 
the particle motion as projected onto the xy plane [3]. 
The in-plane velocity v is constant and the scattering 
probability distribution is oc cosp (Lambert’s law with 



FIG. 1: (Color online) Trajectory of a particle moving bal- 
listically in a circular cylindrical cell with vertical a-axis and 
radius R between successive wall collisions at a and b. A weak 
uniform magnetic Larmor field Bq and a strong uniform elec¬ 
tric field E are applied along, or opposite to, the 2 direction. 
The path segment i shown can be characterized by the wall 
intersection angle p or by the angle where 2ag,i is the 
angle subtended by the chord as seen from the center O of 
the cell. 


in-plane scattering angle p), independent of incident an¬ 
gle Pine and, therefore, independent of the path history. 
Even a single UCN reflection will establish isotropy of 
flight direction. This is also a reasonable model for ^®®Hg 
co-magnetometer atoms which stick to the wall briefly 
before re-emission in a random direction. 

Summarizing previous results |19] we start from the 
Hamiltonian for a classical field, 


Ti = —per ■ B = 


h 

2 


Wg S* 
S —Wg 


( 1 ) 


where (Ta,, ay and Uz are the Pauli matrices, wg = 
—2/rRg/fi, is the Larmor frequency in the vertical static 
magnetic field Bg and p is the magnetic moment of the 
particle. The measurement cell is also exposed to a uni¬ 
form vertical electric field E which gives rise to the mo¬ 
tional magnetic field E x v/c^. 

In Eq. Q, E(<) = wg {Bx{t) iBy[t)) /Bq is the time- 
dependent perturbation due to the small horizontal mag¬ 
netic field [BxiBy) seen by the particle moving from 
point a to point b (Fig.[^ at constant v. In the case with 
small uniform vertical gradient dB^/dz, first analyzed 
in [3], the perturbing field is given as B^ = —Bq^pQ + 
(x/R), By = —(Boy/R, where we use the dimensionless 
parameters Q = v/{Rojq), ( = {R/{2Bq)) dB^/dz and 
p = RujqE/{Bqc^). 

For convenience we define a dimensionless time scale 
T = wgt, reset the clock at the center of each chord and 
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solve the Schrodinger equation for Hamiltonian Q in a 
coordinate system rotating about the vertical z axis at 
the Larmor frequency coq. For Bq pointing up, Wq > 0 
(< 0) for /X < 0 as for neutrons (for /x > 0 as for ^®®Hg). 
Using a second-order perturbation approach we obtain 
the spin-up and spin-down components, ar{T{,,Ta) and 
Pr{Tb, To), of the wave function for a particle starting at 
point a at time Ta = —5i with spin up and measured at 


b at Tb = +6i: 

ar{Si, —Si) = e”*"^ 



/3.(5„-50 = -^( 

Ei„t(5.)-Eint(-50), 

(2) 

where Si = (sinag_i)/U and (see Eq. (9) of [TS]) 


Eint(i) = ^ 


(3) 


For circular cylindrical symmetry, vertical electric and 
magnetic fields E and Bq and uniform vertical gradient 
dB^jdz, the functions v and j3r in ([^ are given by 


Choosing point a (Fig. as the starting point of chord 
1, we set Ta.i = —and obtain for the initial transfer 
matrix Mr^\Si,Ta^i) in ^ 

mP{Su-Si)=Mt{Su-Si) (8) 


with 


Mt((5„ -50 = T{-5,)Mr{5,,-5{)T[-5i). (9) 


The unitary evolution matrix along any chord i, in the 
rotating system, reads 


Mr{6i,-Si) 


/3r(Si,-Si) a*{Si,-6i) 


( 10 ) 


and the transformation matrix from the laboratory to the 
rotating system at point a, 


n-s.) 


- g-z5./2 Q 

0 > 


( 11 ) 


— 3sin5i (?x^ — 5j^) cos 5i-I-2?xi5i sin5i |, 

—5i) = xCU (mi sin5i — St cos 50 = 
xU|77sin5i-I-C (l + sin5^ - 5^ cos 5,j |, (4) 

as shown in [TS], Eq. (42). In Q, ?xi = ?x -f 1, ix = 
( 77 /C) + (cosag,0/ll- 

We write the wave function at the end of n consecutive 
chords in the form 


4^\Sn) = 


(n) 


{Sn) 


= M/")(5„,ra,0 




cos(do/2) 

sin(6»o/2)e*‘^“ 


(5) 


is the same as that for the reverse transition at point h. 

When we perform the matrix multiplications in (§ re¬ 
peatedly for consecutive chords we obtain a sequence of 
matrices m/"^( 5„,—5i) whose general form can be in¬ 
ferred by induction. Just as for a smooth wall the result 
is a unitary transfer matrix of the form 


M/")(5„,-50 


gr -K 

hr g* 


( 12 ) 


for any n> 2. As shown in m. Eq. ( [M| ), the cumulative 
phase advance at the end of chord n, relative to the ro¬ 
tating frame and averaged over initial azimuthal spinor 
angle <l>o relative to initial flight direction, is given solely 
by the argument of gr'. 


{Sipo^n} =-2aiggr zz2lmgr, (13) 


where the initial spinor points in an arbitrary direc¬ 
tion given by angles do (polar) and $0 (azimuthal). In 
the experiments we have do — 7r/2. The recursion rela¬ 
tion linking chains with n — 1 and n segments is (Eq. (33) 
of [H]) 

M(’^)(4,Ta,i) = MT(5„,-5„)i?„M("-i)(5„_i,T,,i), 

( 6 ) 

where 


where the last form in (131 holds since gr differs from 
unity only slightly. 

For any n> 2 we obtain from recursion relation ([^ 


n— 1 n 


X e 


no = l ni=no + l 


(14) 


R 


n 


edas,r.-l-l-ag,„)/2 Q 

0 Q-Kag.n-l+ag.r)/^ 


(7) 


is the rotation matrix, around the z axis, for angular 
change ag^n-i + C(g,n between the coordinate systems of 
segments tx — 1 and n. For each chord the y-axis is chosen 
to lie in the direction of motion. 

In (§ and Q we have left the starting position on 
chord 1, given by arbitrary. The asymptotic fre¬ 
quency shift, for n ^ 1, does not depend on this position. 


where i^i = v{Si,-di) and /3r,i = /3r(5i,-5i). 

In general, reversing the order of two successive chords 
leads to a different phase advance. This property is re¬ 
flected in the form of g* ■ Expression (14) is valid for any 
uq, rii in the range 1 < rxo < txi < tx. 

Using (14) the overall phase shift (13) reads 


{S(po^n) = 2 E ^(^» ~^») + Q (1^) 
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with 


B. Two averaging procedures 


n —1 


Q = - 2 Im| Y (16) 

no = l ni=no + l 


X e 


As shown in [3] for circular cylindrical geometry and 
specular reflection, the frequency shift X(ag,f2) is av¬ 
eraged over all possible chord angles ag using a weight 
factor sin^ ag, thus 


Dividing by the elapsed time To_>„ we obtain 

the average frequency shift {Suj)Q^n- The averaging is 
subtle as described below in Sec. mm after a closer look 
at the definition of angular advance 2Q!g around the cell 
per chord. 

In Eq. Q we have assumed that, going from chord n—1 
to n, the net rotation angle to the new direction of travel 
can be written ag n-i + OLg This is justified if, as in 
[aiis], we define ag&s positive for passage on the right of 
the center O (Fig. Q (corresponding to counterclockwise 
rotation around the cell), as negative for passage on the 
left (for clockwise rotation), and if both chords pass O on 
the same side (which is always true for specular reflection 
and circular geometry). 

However, this definition for ag (which we will later call 
ag,oid) fails if in case of diffuse scattering the consecutive 
chords pass on opposite sides of O. If so the final direc¬ 
tion of travel is opposite that given by ag^n-i+ctg^n- This 
is evident for the special case ag^n = —cig,n-i where the 
particle reverses its direction at the wall reflection point 
and thus has experienced a net rotation by angle tt rather 
than by ag^n^oid T — 0- 

We avoid this inconsistency as follows: For passage on 
the right we define ag as before: for circular geometry, 

ag = p + nj2, for p < 0. (17) 

For passage on the left (p > 0) we choose, instead of the 
negative old value Og^oid = P — 7r/2, continuation of p7| ) 
to the range 7r/2 < Og < tt. Thus 


Og = p -I- 7r/2, also for p > 0. 


(18) 


(181 can also be written ag = Og^oid + tt. With this 
choice of a continuous function for Ug, instead of the 
old function which changes abruptly from -|-7r/2 to —tt 12 
at p = 0, the net rotation angle can be expressed as 
ttg.n-i + OLg,n slso for passage on opposite sides. With 
this choice there is an offset by 27r if the next chord n 
passes on the left, irrespective of the geometry of chord 
n—1. Asa result of this offset, over the course of multiple 
consecutive chords a cumulative offset angle 2m7r may 
appear, with integer m, but this offset is irrelevant both 
for the analysis and the simulations presented below since 
these depend only on the sine and cosine of angles. 

The choice of a continuous function for Oq throughout 
the entire range from 0 to tt, as in ( [T7| , ( [I^ , is also ap¬ 
propriate for the general cell geometries discussed below 
in section cni 


4 

(A(ag,H))^ = - / dag(sin^Q;g)A(ag,D). 
Jo 


(19) 


We use subscript 1 to distinguish this type of averaging 
from a second type introduced below in (221. 


First we show that the averaging scheme (191 can also 


be derived differently and adapted to arbitrary cell ge¬ 
ometry and rough surfaces satisfying Lambert’s law. 

For circular geometry the take-off angle p shown in 
Fig. is complementary to ag: —p = tt (2 — ag, implying 
cos p = sin Q!g. The quantity to be averaged is a fre¬ 
quency, i.e. the phase dp accumulated along the chord 
divided by the traversal time L/v, where L is the chord 
length. Since v = const., the frequency is determined by 
the phase per unit length, dpjL. A longer chord rep¬ 
resents a larger sample for averaging, thus the chord’s 
contribution to the mean frequency is proportional to 
its length L. Including the Lambert factor cosp we ob¬ 
tain the weight factor pi{L,p) ex Lcosp. For circular 
geometry L = 2R cos p, thus pi (p) oc cos^ p = sin^ ag in 


agreement with the weight factor in (19). 


Changing from integration over ag in (19) to integra¬ 
tion over p we have 




^ [ dp {cos^ p)X{p,Vt), 

tVl J-t,/2 


( 20 ) 


where the integral covers the full scattering range 
—TT 12 < p < -|-7r/2. The normalization constant is 
A^i = 7r/2. 

The measured phases are averages over both senses 
of particle circulation around the cell, forward (fw) and 
backward (bw). Eq. (20) represents the part for forward 
motion (in Fig. from atob). The corresponding con¬ 
tribution for backward motion from 6 to a is 




1 

— / dp{cos^ p)X{p-TT,-n), 
JVi J-Tr/2 


( 21 ) 

where we replace p by p — tt and fl by — H to trans¬ 
form to the opposite direction of motion. The fw- 
bw even and odd contributions to the mean values 


are (...)even — (^/^) [(■•■)fw (■■■)bw] and (...)odd — 

(l/2)[(...)fw - (...)bw], respectively. 

Certain averages in ( |15| ) do not involve relative quan¬ 
tities like the frequency (given by the phase per unit 
length) but properties of the chord as a whole, such as 
its length L or angles ag and p or the traversal time 
L/v = 2S/uJo- The contribution of each chord to these 
averages is independent of its length L, thus the weight 
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factor changes to P 2 {p) oc cos p = sin Og. For these quan¬ 
tities the equivalence of averaging procedure (20), (21) 
is, for circular geometry. 


:,odd 


1 

= TT / dp {cos p)Y{p,n), 

tV2 J-Tr/2 

(nP,^)) 2 .^^ = (^(P-^))2.eve„ " ^))2,odd 

1 

= — / dp {cos p)Y{p-Tr,-n), (22) 

^''2 J - Tr /2 


with N 2 = 2. For instance, { 6)2 = { 6 ) 2 ,even = 7r/(4r2). 

In summary, we use averaging 2, given by ( |22[ ) for 
circular geometry, to average properties of the chord 
as a whole, such as <5 or l3r{5,—5) or (as in 

Eq. ([^). We average relative quantities such as fre¬ 
quency or (sin 5 )/6 using method 1 which is given by 
(@, (13) for circular geometry. 

In Sec. Ell we generalize these averaging prescriptions 
to a generic cylindrical cell geometry with arbitrary cross 
section. 


C. Determining the mean frequency shift 

1. For specular wall reflections 

For specular reflection and circular geometry all con¬ 
secutive chords are identical (ofq.i = ag, Si = S, Vi = v, 
/dr,i = Pr) and Eqs. (15), (16) for the net phase shift 
reduce to Eq. (35) of jl9| : 

{Spo^n) = -2arg5r^ 

n—1 

= 2niy{S, —S) — 2s^^^(n — k) sin {2k{ag — i5)) 


k=l 


n>l 


2niy{S, — 6 ) — ns^ cot(ag — S) 


(23) 


with = \l3r{S,—5)\'^. The corresponding asymptotic 
frequency shift (for n 1) is Eq. (37) of [19]: 

{Suj)n-».l ^ / {Spc^n) \ 

Wo \ 2nS / ^ 

V(d,-J)^ _ /P* {5,- 6 )MS,-5) 


(24) 


25 


cot{ag — 5) 


where all mean values are fw-bw averages referring to 


frequency and, therefore, are of type 1 given by (20), 
( [^. The results of [19] , based on (|24| , agree with those 

oTjllTn] 


and given by Eq. (14). We summarize the structure of 
Pr as follows. 

In (14) the functions flr appear in the form of products 
Pr.noPM where uq and ni > ng refer to different chords. 
The following analysis is significantly simplified by noting 
that for completely diffuse scattering the chords are sta¬ 
tistically independent of one another. There is no mem¬ 
ory of the previous path. 


This is an essential point and argued as follows. For 
any trajectory emerging from a given surface element, 
such as at point a in Fig. the 2D Lambert law ensures 
isotropy of take-off direction p, independent of the inci¬ 
dent angle. While this is true for any cell geometry, the 
circular geometry is unique in ensuring a second type 
of uniformity which is also required to render any two 
generations of chords (labeled by indices i and j) inde¬ 
pendent of one another: a circular surface is illuminated 
uniformly by a diffuse point source positioned anywhere 
on the surface. Rays leaving point a at angle p within dp 
intercept a surface strip of width ds = L dp/ cos p where 
L = 2Rcosp is the chord length. Thus ds = 2Rdp is 
independent of impact position. In consequence, at any 
step in any sequence of chords the wall is a uniform source 
of next generation chords. The orientation of the previ¬ 
ous chord does not matter. We will see in Sec. uni that 
this is not true for geometries different from circular and 
in those cases we will have to resort to approximations. 
We obtain an exact solution only for circular geometry. 

In the latter case, the average of a product of chord 
characteristics such as flr, oig, 5 for different chords i and 
j, or of any functions Fi, Gj thereof, can be factorized: 


{F,Gg) = {F/){Gg). (25) 


In the case at hand, where ni ^ no, this applies to the 
product of functions and = 


/3r.nie*^“3'"i in (14). Since all chords are equivalent 


we obtain the same averages for any ng and ni: 


= (/3,e*(“»-^))2, (26) 


where we use averaging scheme 2 since these functions 
refer to a chord as a whole, as defined before (22). 


The same product rule applies to the terms in 
the sum over k in (14) where ng < fc < 
ni. These terms appear as products, such 


P2.i(ag 


H-<5. 


, + i)g 2 i(ag 


rro+2 


-5n 


as 

)+2)02z(ag,TTo+3 —■^TiQ+s) 


average , where we define 


2. For perfectly diffuse wall reflections 

From Eq. the overall phase shift is determined 
by the (1,1) element pr of transfer matrix Mr"^(i5„, —(5i) 


D = {e^^^°‘o-^'>) 2 . (27) 


For circular geometry the explicit form of this mean 
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value is, from (221: 


D = = (l/2)(Aw + ^bw), 

/.7r/2 

Di^= / 

./ —7:12 

^7r/2 

Dbw= / dag(sinag)e2*(-“»-^(-“»’-^» 

Jo 


> -Tr /2 


dp{cosp) e 


2i(ag(p— tt) —5(p—TT, —n)) 


(28) 


where ag(p — tt) = —ag(p) + tt and S(p — Tr,—fl) = 
+S{p, n). In (28) we have separated the terms for fw and 
bw motion and used the symmetry between ranges p < 0 
and p > 0: ag{-p) = -ag{p) + tt, S(-p,fl) = S(p,fl). 
For circular geometry and also for the symmetric geome¬ 
tries analyzed in Sec. |III B] we have D = = £*bw 

For brevity let X = ag - 5, Xq = ag^ng - Sng, = 
= l3r,no, Pi = Prni- Expressed in these 


a 


g,ni ; 


terms the second contribution (16) to the overall phase 


(15) is given by the imaginary part, Q, of 

n—1 n 


s = - 2 J 2 po E 


i[Xo+Xi+2 ,,.1 Xk] 


no = l 
n —1 

= -2 E A 

"0 = 1 


*ei^o 


"1="0+1 
(/3ie*^^)2 


(29) 


X E 

ni=no + l 


2iXr, 


0+^2 X (e^*-^"0+^)2 X ... X (e^*-^"i-02. 


Each average in the last sum of terms is equal to D. 
Therefore this sum is geometric: 


1 + D + D^ + + = 


2 _ £)"-"o 

1-L> 


and Eq. (29) becomes 


^=-2(/3,V^)2(/3,e*^)2 E 

"0 = 1 


I _ £)"-"0 

l-D 


= -2 


1 - 
l-D 


{P*e^^) 2 {Pre ^^)2 
1 - D 


(30) 


where we have used Eq. (26). The second term in the 


square brackets is negligible for n ^ 1 since the mag¬ 
nitude of D is always significantly less than 1. A simi¬ 
lar negligible correction appears if the particle starts on 
chord 1 at a point different from a (Fig. [^. In the EDM 
experiments using UCNs m the starting point is given 
by the beginning of the period of free spin precession in 
the Ramsey scheme and the particle can be at any point 
of the initial chord at this time. 

We derive the ensemble averaged frequency shift (Suj) 
by dividing each chord phase Sip by the passage time 


2d/uJo- Combining (15), (16) and (30) the forward- 


backward averaged frequency shift for perfectly diffuse 
2D wall reflectivity thus becomes 


(Suj) 




n) 


- EmQi, with Qi = (31) 


Wo \Ho^g,^) / l,e^en 

[(/3.e*^)2(/3;e*^/(5)i]_„ + [{Pre^^/SU{P*e^^h 

l-D 


In (31) we use the arithmetic mean value of the two con¬ 


tributions in the numerator of Qi. To show that this 
is justified we employ the following relation between the 
two averaging procedures 1 and 2: For any function F of 
chord variables we have 


(F/S), = {Fh/{Sh. 


(32) 


For circular cell geometry this follows directly from the 
averaging definitions (20) - (22). Since <5 = (cosp)/D, we 
have 


(,5)2 = (I/ 5 V 2 ) 


rTT /2 


l-Tr/2 


dp (cos = Ni/{nN 2 ) (33) 


n 


with normalization constants Ni = 7r/2, N 2 = 2, and 
dp(cos^ p)- 

I-tt /2 


(^/<5)i = ^ / dp(cosV) 

iVi J_^/2 

O ON 

= — / dp{cosp)F=^{F )2 = {F) 2 /{S) 2 , 

iVi J_t ,/2 2Vi 


both for fw and bw motion. This proves (32). Thus, 


the two contributions in the numerator of Qi in (31) are 


equal and it was justified to use their arithmetic mean 
value. 


Relation (32) also holds for the generalization of aver¬ 


aging procedures 1 and 2 to general cell geometry intro¬ 


duced below in (44) - (47) of Sec. Ill A 2 


Using (|32|) the expression for the fw-bw averaged fre- 

((/3;e*^)2(/3re*^)2)even'l 


quency shift (31) becomes 

(5w) { 1 ^) 2 ,even 


Wo (< 5)2 

(^)2,even 


1 


Im ■ 


1 - D 


Im • 


( 5)2 2(5)2 

(/3*e*^)2,fw(^r-e*^)2,fw + (/5re*^)2.bw(/3re*^)2,bw 

1 - D 


(35) 


Eq. (35) can also be understood as: mean frequency shift 


= mean net phase shift divided by mean net time. 

As an example we give the explicit form of (35) for the 
frequency shift cx rp: 


(5w)^ 

Wo 

<52 = 




W 2 

fit 12 


2iV|(5)2 


Im{(l-D)-iQ2}, (36) 


' dp(cos p) sin[5(p, D)] 

-7r/2 

f ' dp(cosp)sin[5(p,D)]e*[-“«(^)-5(^’^)l'^ 

—7r/2 
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We have used the functions v{5,—5) and fir{5,—5) from 
Q and the subscript rf for the terms oc rf, with 
= {r]‘^n‘^/4){2S-sm{26)) and (^'^ 2 ) 2 ,even/(< 5)2 = 
(ri^/2) [l — ri Ji(2/f7)] where Ji is a Bessel function. The 
factor sin 6 in the integrands of (36) is the part of i3r/{i^) 
which is oc rj. The first (second) term in Q 2 describes for¬ 
ward (backward) motion along a pair of chords. 

In Sec. VB[ Eq. ( [^ , we will rederive Eq. (36) in a 
fundamentally different way. 


3. Comparison between specular and diffuse refleetion 


Equation (35) is valid for perfectly diffuse 2D wall re¬ 
flectivity. Comparing it with expression (24) for specular 
reflectivity and noting identity ([M| we see that the first 


term on the right hand side (RHS) of (35), which rep¬ 
resents the single chord contribution, is the same as in 
(p4|), thus independent of surface roughness. 

ensures that the spinor re- 


The second term in 
mains unchanged at each wall reflection. In the limit of 
specular reflection it reduces to the second term in (24) 
as follows: In this case ag is constant along an orbit, thus 
D = The product of two averages in (35) is re¬ 

placed by the average of the single function 
of ag over a random distribution of Og-values. Writing 
Im [ 1/(1 — = Re /[2sin(Q:g — d)] = 

(1/2) cot(og — (5), applying (|32|) and averaging over the 
fw and bw directions we recover (|24|) from 


Prior to comparing numerical results for the frequency 
shifts we extend the constant gradient field, considered 
so far, to include the possibility of a slight misalignment 
between the vertical E field and the Larmor field Bq. 
This type of imperfection of order 10“^ rad or more is 
hardly avoidable in the experimental situation. 


D. Extension to general gradient field 


In Refs, [m [IHl [H] the uniform gradient field has been 
extended to more general macroscopic fields B = Vy de¬ 
rived from a power series expansion of a magnetic po¬ 
tential x{x,y,z), including terms of higher than second 
order, e.g. xy, x'^, x^y, etc. Since the functions v{5,—5) 
and j5r{5,—5) can be determined analytically as in [T^ . 
we can use Eq. (35) to find analytic solutions also for 2D 


diffuse reflection. 

We discuss briefly the effect of tilting the Larmor field 
slightly away from vertical by adding to Bq a small uni¬ 
form horizontal component Bpo = {B^o, Byo). Slight 
misalignment between the volume averaged Larmor field 
and the symmetry axis z is inevitable in the experiments 
and has been taken into account in [iniisa with the result 
that the precession about the tilted z' axis at frequency 

uJq = + B'^q/Bq gives rise to the same frequency 

shift relative to Wq as for a vertical Larmor field with mag¬ 


nitude given by the field magnitude Bl = y Bq + B^q. 

We can derive this result by exactly solving the 
Schrodinger equation for Hamiltonian ([^ in the xyz lab¬ 
oratory frame. The two spinor components a(t), /3{t) are 
determined from m 

a{t) + ^ujQa{t) = 0 , 

R 'I (2*d(t) - uJoait)). {37) 

^o(dJxO ~ t-DyQ) 

The solutions a{t) = cos{oj’t/2) — i(a;o/wo) sin (a;T/2), 
/?(t) = i [(Bxo + *Sj,o)/73l] sin (a;'t/ 2 ), for initial spin 
up, describe uniform precession about the z' axis in the 
xyz frame. In the limit of long time of free precession, 
WQt/(27r) ^ I, which is always realized in the experi¬ 
ments, the precession angles about the z and about the 
z' axis become equal and are given by the field magnitude 
Bl- 

As a result, all results for frequency shifts relative to 
Wo, presented below for purely vertical volume average 
Bq of the field, also hold for a field with slightly tilted 
volume average, except that in this case the shift is mea¬ 
sured relative to Wq. The redefinition of D in terms of Wq 
rather than wq, D —)■ v/{RuJq), implies a negligible scale 
change in all plots of frequency shift vs. D. 


E. Results for circular cell geometry 


To separate the frequency shifts (x (the E-odd shift 
for electric field reversal E —>■ —E, which we normalize in 
the form /{ujo^r])), the second-order motional 

shift cx E^ (denoted by (Suj) ^2 / {ojov'^)) and the second- 
order magnetic shift oc (denoted by (( 5 w)^ 2 /(woC^)) 
we select from Eq. (24) (for specular reflection), and from 
Eq. “ 


(for diffuse reflection) the corresponding terms 
oc Cy, oc 77 ^, and, for E^, those oc 

Figures to 1^ show the results plotted against D. We 
use parameters based on the EDM experiments at the 
ILL [TH3]: E = 10® V/m, Bq = 1 /rT, \{dBo/dz)\ = 1 
nT/m, averaged over a measuring cell with radius R = 
0.235 m and height El = 0.12 m. For these values, |/'| = 
0.00012 and, for i®®Hg co-magnetometer atoms, \r]\ = 
0 . 0001 . 


The integrations required in (35) can be performed an¬ 


alytically, the results involving Bessel and Struve func¬ 
tions. However, most of the calculations were performed 
numerically. 

In Figs. to 1^ the specular data are the same as those 
of pillHj except that in Fig. 2 of m the dependence on fi¬ 
nite n had been shown explicitly and the E^ shift had not 
been averaged over the fw-bw senses of circulation. The 
sharp resonances for specular reflection appear blurred 
since we used a slight imaginary offset of Ug, away from 
the poles, in the numerical integration over ag. 

The limiting behavior for D <C 1 is shown more clearly 
in the right panels of Figs. and where the frequency 

























n 







FIG. 2: (Color online) Upper panels: Normalized frequency shift for if-field reversal (oc E^) as a function of = v/{Rloo) for 
circular cell geometry. The data for perfectly diffuse 2D reflectivity (blue solid line) are compared to specular wall reflection 
(red dashed curve). On the right panel the shift is divided by to show the asymptotic behavior in the adiabatic regime 
II ^ 1. - Lower panels: Comparison of the calculation for diffuse reflection with simulations (orange dashed). 


shifts cx EB and cx have been divided by The 
lower panels of Fig.j^show the agreement, within numer¬ 
ical uncertainty, of i?-odd shift values calculated from 
Eq. with the simulations described below in Sec.lTV 


We observed full agreement also for the second-order 
shifts oc E^ and (x In Fig. the asymptotic form at 
large fl is shown more clearly in the right panels, where 
the functions have been multiplied by 

Comparing the data for smooth walls and rough walls, 
the sharp resonances for specular reflection are largely 
washed out by roughness, as expected and previously 
shown by computer simulations in [5]. 

Both the specular and diffuse reflectivity data in Fig.[^ 
show agreement of the non-adiabatic limit —>■ oo for the 
E-odd shift with the universal Pignol limit [[E], Eq. (5)] 


\ V^o / RB{ 


■{xB^+yBy). (38) 


In circular geometry, the volume average is {xB^ + 
yBy) = -CRBo/2, thus {6oj/{CrjuJo)f^~'' -1- 


III. EXTENSION TO GENERAL CELL SHAPE 
A. Generic cell cross section 

Figure shows a cell of generic cross section. As in 
Fig.[T] we show a trajectory of length L extending from 
point a to point b where a is at position s measured along 
the perimeter from an arbitrary starting point s = 0. The 
chord has the perpendicular distance a from the cross 
sectional center-of-mass point CM, thus a is the variable 
replacing the closest distance x = Rcosag for circular 
geometry, a is positive (negative) for trajectories passing 
the CM on the right (left). 

As in Fig. 2Q;g is the angle subtended by the chord 
L as seen from the CM. It equals the angular advance 
around the cell when the particle travels from a to b. 

For circular geometry we had arbitrarily chosen the ra¬ 
dius i? as a convenient reference length for the definition 
of 11, ( and 77: fl = v/{Rujo), C = RidB^/dz)/(2 Bq) and 
rj = RujqE/{Boc^). For the rectangular cell shape ana- 
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FIG. 3: (Color online) Comparison of normalized second-order motional frequency shift (oc E^) for diffuse reflectivity (shown 
by the blue solid curve) with that for specular reflection (red dashed line) as a function of fl, for circular cell geometry. On the 
right panel the shifts are divided by to show the asymptotic behavior in the adiabatic regime 11 <C 1. 



FIG. 4: (Color online) Comparison of normalized second-order magnetic frequency shift (oc B^) for diffuse reflectivity with 
that for specular reflection as a function of for circular geometry. On the right panel the shifts are multiplied by to show 
the asymptotic behavior in the non-adiabatic regime fl S> 1. 


lyzed below (and also for the equilateral triangle analyzed 
in [20]) it is convenient to choose one half of the average 
side length as the unit of length. Thus, a and L will de¬ 
note normalized lengths which, for circular geometry, are 
given as cr = cos ag = — sin p and L = 2 sin Ug = 2 cos p. 

As before we measure the wall intersection angle p from 
the surface normal and define it as positive for the coun¬ 
terclockwise direction (which corresponds to clockwise 
circulation of the particle around the cell). The trajec¬ 
tory at angle pcM, which passes through the center point, 
divides the region on the right of CM from that on the 
left. We transform from forward to backward direction 
of motion by replacing +n by —H and p hy p — tt. This 
operation reverses the direction of travel from a to b (in 
Figs.[^and[^ to that from b to a, retaining the sign of a 
and 6, but reversing the sign of ag (and adding an offset 


angle tt, see The elapsed time is positive for either 

direction of motion. These symmetries are the same as 
for circular geometry. 

1. The role of chord asymmetry 

As viewed along a chord, a constant gradient field cen¬ 
tered at CM is symmetric about the epicenter c, the point 
closest to CM. Thus, measuring time r from point c we 
have By/B^ = —(V,t and B^/Bq = —(o' = const. 

In general, the times at the chord endpoints, = Sb 
(> 0 for the chord shown) and Tq = —Sa (< 0), will not 
add to zero. As for the circular case we define the average 
as 5 = {6b + 6a) 12.. The imbalance is 6ba = {6b — 6a) /2. 
We also distinguish between two angles subtended from 
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FIG. 5: (Color online) Measurement cell of generic shape sub¬ 
ject to the restriction that any chord intersects the perimeter 
at two points only. The initial point a of chord ah is at loca¬ 
tion s which is measured along the perimeter. As for circular 
geometry, 2ag designates the angle subtended by the chord 
as seen from the areal center of mass CM but the triangle 
a b CM is, in general, asymmetric and the wall intersection 
angle p is no longer directly determined by L and a are 
the length and “lever arm” of the chord. L is divided into two 
sections of, generally, different length: from a to the epicenter 
c and from c to 6 with associated angles ag^a and Qg.t. 


the center: ag^a for the chord segment from a to c, and 
ag^b = ‘2‘OLg — ag^a for that from c to b. 

The imbalance 6ba appears in the expressions for 
v{5b,-5a) and Pr{Sb,-Sa)- 


v{5b,-5a) = —-f 5'^)5 

— 3sini5 — i5^)cosd-|-2Mii5sin(5 | 


/ sin 2(5 . 2 


ha") 


(39) 


(311 and (35) by the replacements 


-^0 — ^riQ — ^riQ 

^ ^1 — 0^5,a,ni ^a,ni 






In the limit \ Eq. ((M|) becomes 


S = —2n 




(42) 


1 - D 

and the fw-bw averaged frequency shift takes the form 


(M 

Wo 


\^/ 2 ,even. 


1 


■ Im 


((/?*e-^°) (/3,e-^0)3 

l-D 


}■ 


(43) 

As in Eq. (351 for circular geometry, the first term in 
(43) is the single chord contribution. The second term 
ensures that the spinor remains unchanged at wall reflec¬ 
tions. The two terms are averaged using the adaptation 
of methods 1 and 2 to general cell geometry described 
in the next section. We will also describe method 3 for 
averaging the product FG in (43), with F = /3*e®^“ and 
G = . For circular geometry this product had been 

averaged as in ( [^ . Its extension to general cell geometry 
is denoted ((F)(G))3. 


2. Averaging procedures for general 2D cell geometry 


For circular cell geometry the surface coordinate s is 
irrelevant and does not appear in the averaging proce¬ 
dures, (20), (21) for method 1 and (22) for method 2. In 


an extension to general geometry we make use of the fact 
that, at equilibrium, all wall reflection points s and chord 
directions are equally probable, as for a 2D gas. There¬ 
fore, at any time the probability, per unit length and unit 
time, for a trajectory to emanate from surface position s 
is the same and the take-off angle p obeys Lambert’s dis¬ 
tribution law (cos p)dp. Thus, averaging procedures (20) 
- (22) are generalized to become (giving only the fw-bw 
even part in each case): 


1 f 

((■■•))l.eve„ = ^ / , 

t/perimeter J—-kI 2 


/ perimeter 
^ [(■•■)«,“1“ ("Osip— TT,— n] : 


dp (cos p)L(s, p) 
(44) 


I 3 r( 5 b, -da) = I sin(5[C (^ + ^) + d] 

— C(5 cos i5-I-sin d|, (40) 


((■•■)), 


1 


/• /•r/2 

/ / 

./perimeter J—-kI 2 


- 7 r /2 


dp (cos p) 


^ [("OsiPin “f ("OsiP—TC, —n] 


(45) 


with ui = 1 -|- (p/C) + (t/D. 

The corresponding change of the transfer matrix ele¬ 
ment g* is taken into account in Eqs. (16), (29), (30), 


with normalization constants 

fr/2 

f perimeter J —t: j 2 


r rr/'2 

Ni= ds dp (cos p)L(s, p) = 2 ttA, (46) 

g perimeter J —t: !2 
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r r '^/2 

N 2 = ds dp (cos p) = 2P 

^perimeter J—7tI2 


(47) 


where A is the cell area and P is its perimeter. 

Expression 7Vi = 27rA holds since, for reversed order of 
integration in (46), the integral f ds (cos p) L(s,p) rep¬ 
resents the totalarea A for all parallel strips of width 
da = ds cos p, extending from the boundary on the left 
to that on the right in Fig.[^ The factor 27r in (46) is the 
result of integration over a uniform distribution of strip 
directions. 

To average the product (/3*e’‘'^°) inEq. (43) we 

have to take account of the fact that consecutive chords 
i and i + 1 are not independent of one another. The 
endpoint of chord i (which is at wall position s) coincides 
with the starting position of chord i + 1 . 

We can quantify this correlation as follows: Consider 
trajectories starting from a given point s in random di¬ 
rection. Except for the singular case of circular geometry 
discussed in Sec. |II the next wall impact point will be 
unevenly distributed over the wall surface and this non¬ 
uniformity extends over a number of consecutive reflec¬ 
tions. This number increases with higher cell asymmetry, 
as for a rectangle with larger aspect ratio. 


As a result, the averaging method used in Eq. (35) will 


be exact only for circular geometry. This is borne out by 
comparison of results for non-circular cells, based on (35), 
with the simulations described below in Sec. It turned 
out that a closer approximation was obtained as follows. 
We average the product of functions E(s,p, O) = /3*e*^° 
and G(s, p, ft) = adopting a third method where 

we treat all pairs of chords as if they were consecutive 
with a common wall collision point s. This model leads 
to 


{F(s, p, n)G{s, p, n ))3 = ^ [ ds 

^■^^3 c/perimeter 

X |[ / dp (cos p)F(s,p, ft) 

f dp (cos p)G(s, p, ft) -I-[fw —>• bw] I (48) 

J—k/2 ^ 


with = 4:P. In Eq. (48) we average over the fw-bw 


directions using the fw—^-bw transformations ft —>■ —ft, 

P~tp— TT, (T— (Ug, 0!g,b) ( — Ug + TT, —Oig^b + 

TT, —ctg,a+'!^), (<5, Sa, Sb) —>■ (S, 6b, 6a) as for a time-reversed 
sequence of chords. 

As always, we obtain the frequency shift by dividing 
the net phase by the elapsed time and obtain Eq. (43) 
via identity Eq. (32). 


Averagin g sch e mes (44) - (47) are adapted to rectangu¬ 
lar cells in ( |A8[ ). For circular geometry, averaging 

procedures 1, 2 from (44), (45) reduce to (20) - (22); 


2((/3*e*'^“) (( 6 re '‘^^))3 reduces to (/3;0*^)2,fw(/3re*'^)2.fw 


(/3*0*^)2, bw(/3re^^2.bw and thus the frequency shift (|43|) 
reduces to Eq. 


II 



FIG. 6: Rectangular measurement cell with side lengths 
2(1 -)- w) and 2(1 — w), thus aspect ratio (1 -I- w)/(l — w). 
Among four groups of paths originating from side IV we show 
a path of group 1 with endpoint on side I (see the text). 


B. Rectangular measurement cell 


Figure shows a rectangular cell with side lengths 
2(1 -I-w) and 2(1 — zc), where w is in the range 0 < zc < 1 
and determines the aspect ratio (1 -|- zz;)/(l — zu). We 
use one eighth of the circumference as the unit of length. 
Referring to the analysis of this system in Appendix 
we use Eqs. (Al) - (|A9|) in approximation (|43|) for the 


frequency shifts and in the simulations of Sec. and 
plot the results in Figs. and for zc = 0 (square cell) 
and w = 0.5 (rectangle with aspect ratio 3). In Fig. 
the values obtained from Eq. ( [4^ for zc = 0 and diffuse 
wall scattering are seen to agree with the simulation data 
within a few percent. 

In Fig. for zc = 0.5, we also compare simulations 
and analysis for specular walls with those for rough walls. 
While the adiabatic and non-adiabatic limits of the sim¬ 
ulations are independent of roughness unnmiiiiiiii], 
the resonances seen in the smooth-wall data are reduced 
by roughness. For the E-odd shift, the diffuse wall 
approximation (43) deviates by 6 % (33%) for zc = 0 
(w = 0.5) from the exact result (38) of [TH] which takes 
the form 


/ 6uj \ e^-e 


"^'-^(l + zu^), (49) 


for arbitrary aspect ratio (1 -I- zz;)/(l — w). 

For specular reflection we can compare the E-odd fre¬ 
quency shift data with a calculation based on an exten¬ 
sion of the analysis of relaxation times by Swank et al. 
[H] to frequency shifts, which is described in Sec. jV C In 
Fig.lH top panel, we observe agreement of this analysis 
with the simulations. 
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FIG. 7: (Color online) On the left: normalized frequency shifts from Eq. (431 for a square measurement cell {w = 0) with rough 
walls as a function of fl. On top: Frequency shift for E-field reversal; center: shift; bottom: shift. On the right panels 

the asymptotic behavior for fl ^ 1 or If 3> 1 is shown more clearly by division by (top and center) or multiplication by 
(bottom). For comparison we include the exact results from simulations. 
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C. Further common results for any cell geometry 

In the adiabatic limit SI ^ 1 we can compare the 
second-order frequency shifts oc E'^ and oc for any 
cell geometry, with the universal expression (18) of [2T]. 
As seen in Figs.[^|^[^and[^we obtain, within statistics. 


where time r = is measured from the epicenter c 
(Fig.[|) and the reference length (given by the radius R 
in case of circular geometry) is set equal to 1. 

Setting Sz(t) = Jzij)/Jxy and denoting the cumula¬ 
tive time elapsed since the start of the Ramsey period by 
'Tcum we obtain from (1^ and (54) the coupled ODE’s 


{6oj)n2 0^1 1 2 o 2 

Wo 2 Wo 


(50) 


for all cell geometries and surface reflection laws inves¬ 
tigated. For circular resp. rectangular resp. triangular 
geometry: ? = 1/4 resp. (1 -|- w^)/3 resp. 1/6, in agree¬ 
ment with Eq. (18) of |21) . 

We also recover the universal adiabatic limit 


{6oj) 


E^-E 

Cv 


2woC?7f2^ 


Ev^ dBz 

1^ Ih 


(51) 


for the E-odd frequency shift. This limit is independent 
of the particle magnetic moment and can be described 
[3 as a purely geometric Berry phase. 

Eurther, we have verified the direct relationship be¬ 
tween the shift cx EB and that oc which is based 
on the general identity (81) derived in Sec. VA3 This 
relationship translates into the identity 


Qrj rf- ’ 

which is clearly seen in Eigs. §§0and§ 


(52) 


C0S($ -f Tcum - A$) 
dr 

- ((a + tjn) sin($ -f Tcum - A$), 

= Sz{t) [(Cct -f 7?fl) C0S($ -f Tcum “ A$) 

-f sin($-f Tcum - A$)]. (55) 


In (55) we take account of the change — A$ of spin 
angle, as seen in the new reference system after each re¬ 
flection. A$ equals the sum of all phase jumps = 

C‘g,b,i-i + Oig,a,i between consecutive chords. 

For each chord, we integrate the coupled first-order 
ode’s (55) numerically, with Sz and $ remaining un¬ 
changed at the impact points. To select the E-odd phase 
shift this process is carried out separately for positive 
r] (corresponding to E pointing up) and for rj replaced by 
—rj (for E down). The difference equals and the E- 
odd frequency shift is obtained by division by Tcum- We 
verified that the asymptotic frequency shifts are indepen¬ 
dent of initial spin angle $o and of starting position at 
any point of the cell area, as expected. 

In Eigs. andwe show also the results for the second 
order shifts oc E^and cx B^. 


IV. SIMULATIONS 

We use the Bloch equation in the form of the two 
coupled ODE’s (44) and (45) of Ref. [3], 

^ = Jz:y [oJz: sin($ -|- wq^) - Wy cos($ -|- Wot)] 

^ [Wa, COs($ -I- Wgt) -|- UJy sin($ -I- Wot)] , (53) 

UjZ tJxy 

where w^, = iVoB^jBo, ujy = uJoBy/Bo, Jz, (J^y) is the 
vertical (horizontal) spin component and $ is the expec¬ 
tation value of its azimuthal angle measured in the sys¬ 
tem rotating about the 2 axis at the Larmor frequency 
Wo. As in [3] we can assume = 0 at the start t = 0 of 
the Ramsey measurement period and, since the pertur¬ 
bation is small, that $ remains small while Jxy remains 
almost exactly constant. Jz oscillates with small ampli¬ 
tude about the equatorial plane. 

As in Ref. [3] we adopt the field model with constant 
vertical gradient. Including the motional field, the net 
magnetic field “seen” along any chord is, for any cell 
geometry: 

^ = -(a-rjn, ^ = -Cy = -(^t, (54) 

-Do -Do 


V. EQUIVALENCE BETWEEN THE 
SCHRODINGER EQUATION AND THE 
CORRELATION FUNCTION APPROACHES 

Note that in this section we mainly use dimensional 
quantities instead of dimensionless parameters like C,, f} 
and Q in order to facilitate reference to some previous 
work. The symbol t will represent the independent vari¬ 
able of the correlation functions and the sign convention 
for the Larmor frequency wq is different. 


A. General proof of equivalence 

As the novel method of calculation used above differs 
from the Redfield theory usually applied in this field we 
now turn to the relation between the two methods. 

In the following we use the direct solutions of the 
Schrodinger equation for confined spin-I/2 particles in 
arbitrary fields in m and m to investigate the condi¬ 
tions necessary for the Redfield theory to be valid. 

We consider the solution of the Schrodinger equation 
for arbitrary time varying fields correct to second order, 
Eqs. (34) and (35) of [37] for the spin-up and spin-down 
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FIG. 8: (Color online) On the left: normalized freqnency shifts for a rectangular measnrement cell with aspect ratio 3 {w = 0.5) 
as a function of If. We compare the simulations and analysis for a smooth wall with those for perfect 2D roughness. On top: 
Frequency shift for f?-field reversal; center: shift; bottom: shift. On the right panels the asymptotic behavior for D ^ 1 

or D 3> 1 is shown by division by (top and center) or mnltiplication by (bottom). In the top panels we have included 
also the FI-odd shift calculated for specular reflectivity on the basis of Swank et al. (2012) and described in Sec. |V C| 
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wave functions gr{t,to) and hr{t,to) in the rotating sys¬ 
tem: 

dr to) 

= 1-^/ / dt"e-“«*"s(t"), (56) 

4 7to Jto 

hr{t,to) = -'- f dt'e-*‘"«‘'s(t')- (57) 

^ Jto 


free spin precession in the Ramsey (free induction decay, 
fid) scheme; thus we are interested in the limit T —> oo. 
We rewrite (621 as 

Suj=—lm dt' dt"e““(*'"*") (E*(t')S(t")), 

7o JQ 


/o ^0 

rT 


(63) 


dr(r-r)e-»"(S*(0)S(r)), (64) 


As above, E(t) = lUx (t) + iiOy (t) is the arbitrary time 
dependent transverse magnetic field. 

Eqs. ( [5^ and (57) are reformulations of Eqs. (17) - 
(19) of |19j with times t, f, t" no longer restricted to a 
single chord but considered a common variable. We have 
adapted the terminology of m to that used in [TH] and 
in the present article. In symbols ar and j3r were 
used for the wave functions at any time, not only for a 
single chord. We replace them by gr and as in (12). 


1. Frequency shifts 


where the second step is a well known result [28l [29] . 

We have made the usual assumption that the field 
fluctuations have zero time average and are stationary, 
i.e. (E*(t')S(t")) = {E*{0)E{t = t' -t")) = R{t) = 
R{—t). This also applies to the case of specular reflec¬ 
tions when we take the ensemble average over different 
trajectories (different Og’s) and different starting posi¬ 
tions as has been shown in m- 

In the usual physical situation we have for the obser¬ 
vation time T: T ^ Treiax >> Tcorr where Treiax is the 
spin relaxation time and R (r) —^ 0 for r > Tcorr so that 
in this limit 


The wave function of the spin system, tp (t) is given 
by the matrix (12) acting on the initial state. We take 


that initial state to be the state corresponding to the 
spin pointing in the x direction and evaluate ((T+)^ = 
{ax + in the state tp (t). (cr+)^ represents the rotat¬ 
ing component of the spin as seen in the rotating coor¬ 
dinate system. The time dependence of its phase gives 
the frequency shift and the decay of its magnitude gives 
T2. Thus 


{<^+)r = ( 9 *^ - 


(58) 


{fJ+}r = 1 - 


= 1 -I- 2:2 


(t') / (t") 


(59) 

(60) 


neglecting rapidly varying (2a;o) term. 

We then evaluate the ensemble-averaged frequency 
shift, using 

{S(p) = arg {a+)^ = arg (1 -k Z 2 ) ~ Im [za] 

= --Im/ dt'e“«‘'E*(t') / (t") 


(61) 


{Suj) 


{H)t 

T 


(62) 




dt" 



where we have set tg = 0- T is the total elapsed time 
which, in the experiments, is given by the long period of 


Suj=llm [ dTe^‘^°^ (S*(0) E(t)) 

2 Jo 

1 r°° 

= - / dr coswoT Im (E*(0) E(t)) 

2 Jo 

1 r°° 

+ z dr sin wgr Re (E* (0) E(t)) (65) 

2 Jo 

which is the Redfield result [13]. 


2. Relaxation rate I/T 2 

We now calculate the relaxation rate I/T2. 

Since l^al << 1 

|(a+)j" = l + 2Rez2 (66) 

l(f^+)rl = 1 + Re2;2 

= 1 - i Re / dt' [ dt"e*‘^“(*'"‘") (E* {t") E {t')) 

2 Jto Jto 

(67) 

We now specialize to the case of a stationary stochas¬ 
tic system where (E* (d) E {t")) is a function of {t' — t") 
only. Consider a square region of the t' plane between 
[to^to), {to,t), {t,to) and {t,t). (See (jlH]) for a discus¬ 
sion of this argument). Then the double integral over 
the top half {t' > t") is seen to be the complex conjugate 
of the integral over the bottom half {t' < t"), so the last 
term is given by the integral over the entire square. 

As a result of this we have (again putting t' — t" = r) 
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l(^+).l 


l(^+).l 


1 - i Re / dt' [ (E* {t') E {t")) 

2 Jo Jo 

1 - i Re f dr {t- |t|) (E* {t') E {t' - r)) 

2 7o 

1 - Re f (E* (t') E (i' - r)) (68) 

2 Jo 

1 - (69) 

^2 


where the result is valid for times t > Tcorr- Substituting 


E {t) = ujx (t) + iujy (t) (70) 


we obtain the usual Redfield result: 


1 


1 

2 



(E* (t') E (t' 


t)) 



I 

V 


cos Wo'T { 
+ sin woT 


UJ^ (0 uJx {t-r) \ 

+ UJy it) UJy (t-T) / 

UJy (<) UJx (t-r) \ 

, (0 UJy {t-T) / 


(71) 
\ 

/ 

(72) 


The second term is absent in the usual treatments as it 
is normally assumed that the cross correlation between 
the components of the fluctuating field vanishes but this 
does not hold in the presence of the E x v held. 

Thus we have shown that the Redheld expressions for 
the frequency shift and T2 relaxation rate, see also [14], 
follow directly from the second order solution to the 
Schrodinger equation for arbitrary time dependent held 
with zero time average. 

To calculate the longitudinal relaxation time, Ti, we 
calculate tp (t) starting in the spin up state and and then 
evaluate (cr^),, = \gr\^ — \hr\^ with the well known result 
1/Ti = 2 IT 2 where I/T2 is the transverse relaxation rate 
in the absence of huctuations in i?z, m- 


3 . Expressions for the frequency shifts 


Using (70) we hnd 


E = {-fBx - bvy) + i {-yBy + bv^) 


(73) 


where, as above, we are considering a magnetic held B 
in the presence of a motional magnetic held (E X v)/c2 
produced by an electric held in the z direction parallel to 
B (we use b = and 7 = ojqIBq is the gyromagnetic 

ratio of the particle) 


Re (E*(0) E(t)) = (72 {Bx (0) R, (r) + By (0) By (r)) 

+ b^ (vx (0) Vx (r) + Vy (0) Vy (r)) ), 

Im (E*(0) E(t)) = -bj(vx (0) Bx (r) - Bx (0) Vx (r) 

+ Vy (0) By (r) - By (0) Vy (t) ) 

= 2lrf{Bx{0)Vx{r) + By{0)VyiT)) (74) 


where we have assumed that the motions in the x and 
y directions are uncorrelated. Compare [18], equations 
(8) - (10). Specializing to the case of a linear magnetic 
gradient, dBx/dx = dBy/dy = G, we have 

Re (E*(0) E(r)) = 'y‘^G‘^{x (0) x (t) + y (0) y (r) ) 

+ {lElc^y {vx (0) Vx (t) + Vy (0) Vy (r) ), 

Im(E*(0)E(r)) 

= 2'y'^{GElc^){x (0) v^ (r) + y (0) Vy (r) ). (75) 

Thus according to equations (65) and ( |7^ , see also 
[10] . the frequency shift ~ E^ is given by 

5we^ = (7E/c^)^Im [5™ (wq)] , (76) 

poo 

Sxv{oJo)= dTe*‘"«"5™(r) 

Jo 

with Syy{T) = {vx (0) Vx (r) + Vy (0) Vy ( t )) , 
while the linear in E shift is given by 

Sujeb = 2^^{GE/c^)Re SxviuJo) 

nOO 

with Sxv{uJo) = / Sxv, (77) 

Jo 

Sxv (r) = {x (0) (r) + y (0) Vy (r) ) = f Syy (fi) 

Jo 


poo pT 

thus Sxv (wo) = / dTE‘^°'" / Svv (m) 

Jo Jo 


dy 


1 ^ 


= — / dTE'^°'"Svv (t) = —Svv (wo) 

UJo Jo UJo 


(78) 


z — 1 — 

and similarly, Sxx (wq) = - — Sxv (wq) = -^Sw (wo) ■ 


Wo 


Wn 


Substituting (781 in (771 we obtain 


5loeb = - ^5 Im (wo)l 

and comparing (80) with (76) we hnd 


duJEB _ ^ Gc^ 

Slve^ Eujq ’ 


(79) 


(80) 


(81) 


i. e. the frequency dependence of the two effects differs 
only by the factor ujq. This was noted by Pendlebury et 
al. [3| for the case of a circular cylinder with specular 
walls but our result is universal. It holds for all geome¬ 
tries, held conhgurations and collisional regimes and is 
clearly shown in Figs. iilzl and[^ Using ([78|we can 
also show that the result given in Eq. (64) of 35] does 
indeed agree with previously published results contrary 
to the doubt expressed in [32] . 


B. Correlation functions for a circular cylinder 
with diffuse wall collisions 

We calculate the velocity correlation function for a cir¬ 
cular vessel with diffuse rehecting walls as follows. 
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We consider a trajectory that starts on a chord with 
agfi (Fig. at a distance x from the end of the chord. 
The particle reaches the end of the chord at time t„=o- It 
then travels on a chord ag^i until a time ti. Between 
and tn+i the chord angle is ag^n+i. At each wall collision 
at tn the velocity vector rotates through an angle SO — 
ttg.n+i + o.g,n- For negative Pn+i (motion in the positive 
CCW direction), 69 = 7r+p„+p„_|_i for both positive and 
negative while for positive Pn+i (motion in the CW 
direction) we have SO = —ir + pn+Pn+i- Since tt = —tt so 
long as angles are concerned all that matters is whether 
we have an odd or even number of tt's. This will result 
in a varying sign when we take the sine or cosine of the 
angle. We have that between and tn+i the velocity 
makes an angle On with its original direction with 


COs6'n = (-l)''^^COS (po+/>n+l + 2 ^ pfc ) • (82) 




The wall collisions occur at times 


2R 

tn = - 2^ cos Pk 

^ \fc=l / 


to 


(83) 


with to = x/v. 

So the correlation function, after t = tQ, consists of a 
series of pulses between times tn and tn+i with heights 
(normalized to v^) cos On- 


The Fourier transform Syy (ui) (76) will be a sum 


of terms each corresponding to one of the rectangular 
pulses. Each term is given by the Fourier transform of a 
rectangular pulse. 


Sn (oj) = 2 COS 0, 


sin^^ 

, C '' ' 

LU 

sinSrt 


^ ^ ^...Sn+l -ic. ( *"+(+*" ) 

= 2cos6^n-e V 2 y 


LJ 


(84) 


^ 2 ^-iu^x/v cos On e-»(^..+i+2 s,) ^ 


cu 


where Sn = cos pn and 


{St)n = tn+1 -tn= - (cOS Pn+l) ■ (85) 

V 


We have used, from (83), 

^n+1 


cospn+1 +2^cos/?fc +to- (86) 


k^l 


We have to average the result over all possible starting 
positions and directions of motion. That means we av¬ 
erage over all possible initial chord angles, OLg^, and all 
positions along the chord given by the distance x from 
the end. x varies between x = 0 and x = 2R cos p. Aver¬ 


aging (84) over x gives 

p2R cos pq 


2R cos po Jo 


^ — ILOxfl 


dx 


u) R cos Po 

r- C- . 


- -sin fi? (cos Po) -1 

OS iOn V V J 


By (a.) = 

do w 

We write cos On = (e*®" -I- /2 and consider only 

the first term. We will then add in the second term at 
the end. 

Summing (w) from (87) over all chords except the 
first (labeled 00) we have, using (82), 

_ 2sind0g,(p„_5„) 
uj So 

N 


(87) 


sin S„ 


( 88 ) 


X (-l)"sin^„+ie 


i(p„+i-5„+i) 


n—0 


For j = l...n the angles pj are independent of po, Pn+i 
and of each other, so the average of (881 can be written 


LU 

N 


2 /®^e*(Po-^o)^^sm5„+ie 


i{Py. + l-Sy+l) 


n—0 

2 / sin i5i 

UJ 


2i(pj-5j) 


So 

1 


Ogdpo-io)^ (Jsm5„+ie 


i{Pn + l-Sy+l) 


X 


1 -)- ^e^dpj 

To this we add the contribution of the second term, , 
in the expansion of cos and the contribution of the first 
chord which adds a term equal to for a time from t = 0 
to t = to = x/v and will have a spectrum 

Boo =2 bJ g-»c.(^)_ 


(89) 


Its average over x is 


1 


n2R COS Po 


Soo{uj) — . 

cart cos po 7o 


sin 




(90) 


and the complete expression for the spectrum is 


Syv (^) - 


2 / sin Si 


UJ \ So 

X (sin(5„+ie*^^'‘+^“‘^"+^ 


d gdPo-i5o) 


1 -)- ‘^j)^ 


+ KPj -Sj)^ {Pj - TT - (5j )] 


-k 


1 /I 1 


ujSo 


2 2® 


- iSf 


(91) 
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The average has to be taken over all possible starting po¬ 
sitions and directions of velocity (x, po) ■ Since the num¬ 
ber of particles located on a given chord is proportional to 
the length of the chord (Ij = 2Rcospj) we use averaging 
method 1 for the first chord. However because following 
each wall collision the angles pj are distributed by a co¬ 
sine law (Lambert’s law in 2 dimensions) we use method 
2 for all other chords. Then using Eq. (76) we can write 




{"fvE/c^ 


= Im [Syy (wo)] 


^^sin Sq€ 




= Im ^ liiQ (i5o>2 

+ [(“SJ - ^j) ^ {-<^9,3 + TT - (5j)] 

i (sin2(5o) - (5o 


Wo (^ 0)2 \ V 2 


(92) 


in agreement with the Eq. (36), taking into account a 


different sign convention for wq. We have included the 
normalization factor in (92) and, on the last term, 
used (32) and again (76). 


C. Correlation functions and frequency shifts for a 
rectangular cell with specular walls 


Substituting expression (75) for Im (E*(0) E(r)) in 
equation (65), noting that 

{x (0) (r)) = (0) a; (r)) (93) 


and integrating by parts we find for the linear in E fre¬ 
quency shift [50] : 


Sujeb - Y2 


Ll) 


(94) 


-Wo/ dr sin uJoT{x{0)x{T) + y{0)y{T)) 

Jo 

a result first obtained in [15]. Swank et al. [15] (Eq. (17)), 
using the form of transport propagator derived by Ma- 
soliver et al. [31], have given an analytic result for the 
Fourier transform, S^x (w), of the correlation function 
{x (0) X (r)) in the case of a rectangular cell with spec¬ 
ularly reflecting walls and gas collisions, which together 
with (]^ yields: 


(95) 


Sujeb — + ^1) + 


Im • 




odd 


[{l+iuJoTc)^-\-{n^7rl^y 

+ {x^y) 


P72- 


where lx,y = vTc/Lx^y = \clLx,y, Lx^y are the lengths of 
the sides, v is the velocity, Tc is the average time between 
gas collisions and Ac is the mean free path. This is valid 
for all values of the mean free path, so taking the ballistic 
limit Ix^y S> 1 we obtain the result plotted in Fig. I^with 
reference to Swank et al. [E]. More details can be found 


VI. SUMMARY AND CONCLUSIONS 

Solving the Schrodinger equation directly we have ana¬ 
lyzed the frequency shifts in searches for an electric dipole 
moment of the neutron using ultracold neutrons and co¬ 
magnetometer atoms confined in measurement cells. The 
cells are assumed to have a perfectly flat top and bottom 
and can have rough cylinder walls (Sec. 0 and a cross 
section of generic shape (Sec. Ill A). For perfectly diffuse 
2D roughness scattering on the walls and circular cell 
cross section we derived expression (35) which is valid 


for any in-plane particle velocity v. Its generalization 
to arbitrary cell geometry in the form of Eq. (43) is an 
approximation. 

We have specifically i nvesti gated the geometry of a 
general rectangle in Sec. IIIBj ) and compared all results 
with the simulations described in Sec. EYi So far, general 
results valid in a broad range of particle velocity had been 
reported only for circular and rectangular geometry and 
specular reflection [SKnidSKis]. 

In Sec. IID we showed that a Larmor field slightly 


tilted away from the vertical symmetry axis of the cell, 
due to misalignment inevitable in the experiments, can 
be taken into account as in Ref. m- The shifts remain 
unchanged but refer to the Larmor frequency ojq given 
by the field magnitude, rather than to wq. 

The adiabatic limit (D ^ I) of the E-odd shift agrees 
with the purely geometric Berry phase while those for 
the second-order frequency shifts, oc E^ and oc agree 
with the general results of [5T]. This consistency and, 
generally, the equivalence between previous approaches 
and our analysis are elaborated in Sec. jVj We show by 
direct calculation the agreement of the results from the 
Schrodinger equation with the Redfield theory for the 
cases of a circular cell with diffuse reflecting walls and of 
a rectangular cell with specular walls. 

These comparisons serve to highlight the conditions of 
validity of the Redfield theory which makes use of a num¬ 
ber of assumptions m and has resulted in a considerable 
literature concerning its range of validity, e. g. [22] . Our 
results can be applied to cases where the Redfield re¬ 
sult no longer holds, such as non-stochastic or partially 
stochastic systems. An example of the latter is the in¬ 
evitable situation where the static fields E and Bq are 
not exactly aligned or where, for non-spherical cells, Bg 
is not exactly aligned with the symmetry axis. 

Our method also applies to cases where the observa¬ 
tion time is on the order of or shorter than the correlation 
time, and thus we can describe the transient spin dynam¬ 
ics, i.e. the gradual development of the shift with increas¬ 
ing time subsequent to the start of the free precession in 
the Ramsey scheme. This information is more detailed 
than the transverse relaxation time T 2 derived from the 
Redfield theory or directly from the Schrodinger equation 
[27] or the Bloch equation [32]. Similarly, studying the 
vertical spin oscillations as in [19] can provide informa¬ 
tion not contained in Ti. Our analysis in the present arti¬ 
cle does not include gravitational depolarization [55] |53] 
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due to the effect of gravity on the UCN spectra in a spe¬ 
cific apparatus. 
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Appendix A: Analysis of rectangnlar cell geometry 

To cover all possible trajectories for the rectangular 
geometry shown in Fig. [^it suffices to analyze the four 
groups of chords listed in Table |l] These start from side 
IV (or equivalently, from II) in any direction (—7r/2 < 
p < 7r/2). The paths of group I (4) end on side I (III). 
We have divided those ending on side II into groups 2 
and 3. They are separated by the angle pcm{x,w) = 
arctan [x/{l — w)] for the path passing through CM. An¬ 
gle pi{x,w) = — arctan [(1-I-w — a;)/(2(l — w)] is di¬ 
rected towards the upper right corner. For trajectories 
originating from sides I or III we replace w by —w. 


TABLE I: p-ranges for the four path groups for a rectangle 


Group j 

p-range 

1 

-f < P < pi{x,w) 

2 

pi{x,w) < P < pcm{x,w) 

3 

pcyi{x,w) < p < -p\{-x,w) 

4 

-pi{-x,w) < p < f 


The characteristics for each path are as follows: 

Li{x, p,w) = —{1 + w — x)/ sin p (Al) 

Lj {x,p,w) = 2{1 — w)/ cos p, j = 2,3 (A2) 

L4(a;,p,t/;) = (1-I-ui-I-a;)/sinp. (A3) 


For j = 1,..., 4 we have 


aj{x,p,w) 

cos{2agj) 


(x cos p — (1 — w) sin p) sign(r2) 


(A4) 






sin(2agj) = 


.jvfTw]' 


Vj {x, p, w) = [x sin p + (1 — w) cos p\Lj 
- x"^ - {1- wf, 

Wj{x, p, w) = [xcosp — (1 — w) sin p]Lj, (A5) 


ag^j{x,p,w,V,) = -sign(H) arccos(cos2agj), (p < pcm) 
= ^sign(H) [27r - arccos(cos 2agj)], (p > Pcm); (A6) 

o:g,a,j{x,p,w,n) = (- -h p - pcM)sign(H), 

— ag^aj, 

6j{x, p, w, H) = Lj{x, p, w)/{2Q), 

Sa,j{x,p,w, H) = ^ [(1 - w) cos p -I- X sin p], 

^b,j ~ ^a,j 5 

Sba,j = Sj - Sa,j. (A7) 


The signs appropriate to the different ranges of p and to 
backward vs. forward motion are shown explicitly. For 
the arccosine function we use the range from 0 to tt. 


Adapting Eqs. (441 to (47) to rectangular geometry we 
average the frequency shift (43) over the four groups as 




pl+lil 


dx i 


dp cos p 


^— 1—w ^ p—range 


^ [("Ofw (--Obw] “1“ ^ 


(AS) 


with normalization constant N 2 = 4, and the product 
term in (43) in the form 


1 


4 i<l+w 


dx / dp cos p 

— 1 —Lu range 


^ ([-^Ifw [G'lfw + [-^Ibw [G]bw) + («^ -w) |, (A9) 


with N 3 = 8. Some averages can be determined analyti¬ 
cally, e. g. { 6)2 = 7r(l — but most integrations 

had to be performed numerically. 
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